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Abstract 



We perform the baseline and energy optimization of the Neutrino Factory including the 
latest simulation results on the magnetized iron detector (MIND). We also consider the 
impact of r decays, generated by z/^ — )■ Ur or i/g — )■ appearance, on the mass hierarchy, 
CP violation, and ^13 discovery reaches, which we find to be negligible for the considered 
detector. For the baseline-energy optimization for small sin^2^i3, we qualitatively recover 
the results with earlier simulations of the MIND detector. We find optimal baselines of 
about 2 500 km to 5 000 km for the CP violation measurement, where now values of as 
low as about 12 GeV may be possible. However, for large sin^ 26'i3, we demonstrate that the 
lower threshold and the backgrounds reconstructed at lower energies allow in fact for muon 
energies as low as 5 GeV at considerably shorter baselines, such as FNAL-Homestake. This 
implies that with the latest MIND analysis, low- and high-energy versions of the Neutrino 
Factory are just two different versions of the same experiment optimized for different parts 
of the parameter space. Apart from a green-field study of the updated detector performance, 
we discuss specific implementations for the two-baseline Neutrino Factory, where the con- 
sidered detector sites are taken to be currently discussed underground laboratories. We find 
that reasonable setups can be found for the Neutrino Factory source in Asia, Europe, and 
North America, and that a triangular-shaped storage ring is possible in all cases based on 
geometrical arguments only. 
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1 Introduction 



Neutrino oscillation experiments have provided compelling evidence that the active neutri- 
nos are massive particles pointing towards physics beyond the Standard Model. In a 
three-generation scenario, there are two characteristic mass squared splittings {Am'^^ , Amgi) 
and three mixing angles {9i2 , 6^13 , 623) as well as a CP violating phase 6cp affecting neutrino 
oscillations. A global fit of solar, atmospheric, accelerator and reactor neutrino oscilla- 
tion experiments, yields the following parameter ranges for the oscillation parameters at 
the la level [2]: Am^^i = (7-59 ± 0.20) x 10"^ eV^ jAm^J = (2.46 ± 0.12) x lO'^ eV^, 
012 = (34.4 ± 1.0)°, 623 = (42.81^:^)° and ^13 = (5.6^|?)°. Even, within the three flavor 
framework, there are still unknowns: the mass hierarchy (MH)-Am|]^ > (normal ordering) 
or Am|^ < (inverted ordering); the value of and whether there is CP violation (CPV) 
in the lepton sector. The experiment class, which may ultimately address these questions, 
is the Neutrino Factory j4|[5]. 

In a Neutrino Factory, neutrinos are produced from muon decays in straight sections of a 
muon storage ring. The feasibility and a possible design of a Neutrino Factory have been 
subject of several, extensive international studies, such as in Refs. [6]-[8]. The International 
Neutrino Factory and Superbeam Scoping Study (sjjlO] has laid the foundations for the 
currently ongoing Design Study for the Neutrino Factory (IDS-NF) [5]. The goal of the IDS- 
NF is to present a conceptual design report, a schedule, a cost estimate, and a risk assessment 
for a Neutrino Factory facility by 2013. The IDS-NF defines a first-version baseline setup of 
a high energy neutrino factory with E^^ = 25 GeV and two baselines Li ~ 3 000 — 5 000 km 
and L2 — 7500 km (the "magic" baseline [Tl]) served by two racetrack-shaped storage 



rings, with a muon energy of 25 GeV (for optimization questions, see, e.g., Refs. 12 19]). 
This setup has been demonstrated to have excellent sin^ 26'i3 reaches for addressing the 

to be robust against many potential new 
and to be useful for degeneracy resolution 



open questions in the three flavor scenario 17 



physics effects 19,20 or systematical errors 21 



independently of the finally achieved luminosity 11 ; for the physics case for the very long 
baseline, see also Ref. llSl. 



The appearance signal in a neutrino factory consists of so called wrong-sign muons [e.g., 
from Uf. — )■ Ufj) and therefore a detector which is capable of measuring the charge sign of 
muons is required in order to distinguish this signal from the right-sign {e.g., from z/^ — )■ z/^) 
muon background. The most straightforward solution towards a high fidelity muon charge 
measurement is a magnetized iron detector (MIND). With a MIND, the achievable levels 
of muon charge identification allow for CP violation measurements in the muon neutrino 



appearance channels 13 , 14 



The optimal MIND detector has backgrounds (such as from neutral currents or charge mis- 
identification) at the level of about 10~^ to 10~^, and the potential to measure the muon 
charges at relatively low energies down to a few GeV. The importance of the precise location 



of the detection threshold was discussed in detail in Refs. 17,22 . As the design of the 



Neutrino Factory matures, more refined detector simulations have become available 23 , 24 
especially in comparison to the IDS-NF baseline 1.0 [5]. Compared to the older analyzes 

^ Note, that recent hints for for ^13 > |3j are inconclusive and need to await more experimental data. 
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these new simulations provide the detector response in terms of migration matrices mapping 
the incident to the reconstructed neutrino energy for all individual signal and background 
channels. An optimization of the cuts has lead to a lower threshold and higher signal 
efficiencies than in previous versions, while the background level has been maintained in the 
most recent analysis [24]. In addition, separate detector response functions for neutrinos 
and anti-neutrinos are available, and it turns out that the detection efficiency is better 
than the z/^ detection efficiency, which partially compensates for the different cross sectionsj^ 



The MIND detector has been studied in Ref. 23 , 24 as generic neutrino factory detector 



and a specific detector of similar type is proposed for the India-based Neutrino Observatory 



(INO) to measure atmospheric neutrinos 25 . The detector at INO may serve as a Neutrino 
Factory far detector at a later stage. 



Most recently, the background from r decays was discussed for disappearance 26 and 
appearance [27] channels. These taus arise from charged current interaction of u^- which are 
due to oscillation, e.g., for /x"*" stored: 

App.: z/g — )■ z/t- — 7- ""-T" (background) versus ^ ^ fJ,~ (signal) (1) 

Disapp.: z/^ — )• z/^ -> /i+ (background) versus z/^ — > z^^^ — > /i^ (signal) (2) 

The reason for these muons to contribute to the background is that the MIND cannot 
resolve the second vertex from the r decay, in contrast to OPERA-like emulsion cloud 
chamber (ECC) pSl. In principle, the muons from r decays carry information which may 



be used for the standard oscillation 29 , 30 or new physics |31j measurements 



An alternative version of the Neutrino Factory with respect to the IDS-NF baseline has 



been proposed in Refs. 32 37 . The key difference is to replace the MIND with a magne- 
tized totally active scintillator detector (TASD). The TASD, being fully active, has a lower 
threshold and better energy resolution. The better detector performance and an optimiza- 
tion of the front-end increasing the intensity have allowed a version of the Neutrino Factory 
with i?^ ~ 5 GeV and a baseline possibly as short as L ~ 1 300 km, corresponding to FNAL- 
Homestake. This version is usually called "low energy Neutrino Factory" (LENF) and it 
is found that the LENF has especially good performances for large sin^ 2^13. In addition, 
the performance of a TASD allows to exploit the "platinum channel" (z/^ z/g), however it 
turns out that it is of little practical value j36] . 

The recent simulation results for the MIND have made the performance margin between 
TASD and MIND considerably smaller and therefore, we will show that the distinction 
between the low- and high-energy Neutrino Factory is somewhat artificial and merely cor- 
responds to two extreme corners of a common parameter space. 

The phenomenological discussion of the Neutrino Factory has so far been performed mostly 
in an abstract baseline-energy space. While the energy is a continuous variable, it is not 
obvious that all baselines can be realized from any accelerator site. Therefore, we will present 
a comparison of physics performances for a judicious choice of accelerator and detector 
locations. It seems unlikely that a machine of the size and complexity of the accelerator 
part of a Neutrino Factory would be built on a green-field site and therefore, we assume 



^The difference in neutrino and anti-neutrino response is due to the different y-distributions (24j. 
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that it will be co-located with an existing, large accelerator facility. To be specific we 
consider: CERN, the Fermi National Accelerator Laboratory (FNAL), the Japan Proton 
Accelerator Research Complex (J-PARC), and the Rutherford Appleton Laboratory (RAL) 
as potential accelerator sites [lO]. For the choice of potential detector sites the issue is less 
clear-cut, since, at the current stage, there is little information on the required amount 
of rock overburden for a MIND (or TASD) to perform satisfactorily. Therefore, we make 
the conservative choice and assume that a Neutrino Factory far detector requires a similar 
amount of rock overburden as other neutrino experiments do. Under this assumption, a 
natural choice of candidate detector sites is given by candidate detector sites for other 
neutrino experiments. Fortunately, lists of candidate sites for general neutrino experiments 
have been compiled for the US in response to the National Science Foundation (NSF) call for 
proposals for a Deep Underground Science and Engineering Laboratory (DUSEL) [38] and 
for Europe in the context of the Large Apparatus studying Grand Unification and Neutrino 
Astrophysics (LAGUNA) j39] study. In North America, we consider eight locations: Soudan, 
WIPP, Homestake, SNOLAB, Henderson, Icicle Creek, San Jacinto, and Kimballton. In 
Europe, under LAGUNA, there are seven possible candidate sites: Pyhasalmi in Finland, 
Slanic in Romania, Boulby in UK, Canfranc in Spain, Frejus in France, SUNLAB in Poland, 
and Umbria in Italy. Along with these seven sites, we also consider Gran Sasso National 
Laboratory (LNGS) in Italy and Gran Canaria in Spain. We will complement these lists 
of detector sites by the Asian facilities: the Kamioka mine in Japan, the proposed Chinese 
underground laboratory at CJPL, YangYang in Korea, as well as INO in India. 

This paper is organized as follows: We describe our methods and implementation in Sec. [2] 
After that, we update the simultaneous optimization of baseline and muon energy in Sec.[3]in 
a green-field scenario. In Sec. |4| we discuss the selection of specific sites, the site geometry, 
and the possibility to use a triangular-shaped storage ring. Furthermore, in Sec. [5} we 
quantify site-specific performance of the Neutrino Factory. Finally, we summarize and draw 
our conclusions in Sec. |6j Details for the assumptions for the individual accelerator and 
detector sites can be found in Appendix |Xj The sensitivity curves for all possible considered 
site combinations are given Appendix [B} The individual data files for the curves are available 
for download at Ref. |40] . 



2 Simulation method and performance 

In this section we describe our simulation method and we show the difference to the IDS-NF 
1.0 in terms of event rates. We also compare the performance resulting from the different 
detector simulations, and we compare the performance between one and two baselines. 



2.1 Simulation method 



For the simulation of the Neutrino Factory, we use the GLoBES software 
description of the experiment is based on Refs. 



17,22 



where we use the parameters from 
the IDS-NF baseline setup 1.0 (IDS-NF 1.0) described in Ref. IS) (note number IDS-NF- 



41 42 . The 



002). The detector description of this setup is based on Ref. [9], which has been updated 
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in Refs. 23 , 24 . In this section, we compare these three detector descriptions, whereas we 
use only the most recent version, Ref. [24], in the following sections. IDS-NF 1.0 uses two 
magnetized iron calorimeters (fiducial mass 50 kt) at L = 4000 km and L = 7500 km. There 
are two racetrack-shaped storage rings pointing towards these detectors, with a luminosity 
of 2.5 X 10^° useful muon decays per polarity, decay straight, and year, i.e., 10^^ useful 
muon decay per year. We assume a running time of 10 years, i.e., 10^^ useful muon decay in 
total. The parent muon energy is assumed to be = 25 GeV. The considered oscillation 
channels are: 

appearance: z/g — ^ for /i^ stored , (3) 
appearance: z/g — ?■ t'^ for fi~ stored , (4) 
z/^ disappearance: z/^ — t- z/^ for /i^ stored , (5) 
z/^ disappearance: z/^ — )■ z/^ for /i^ stored . (6) 

Since the luminosity changes if one or two storage rings are required, i.e., one or two 
baselines are operated, and the efficiency of a triangular-shaped ring, which will discuss 
later, is different, it is convenient to re-parameterize luminosity in terms of a scale factor 
(SF) [37]: SF=1 corresponds to the above mentioned parameters 2.5 x 10^" useful muon 
decays per polarity, decay straight, and year. If only one baseline is needed, then all muons 
can be injected in the same storage ring, and SF=2. If, on the other hand, a storage ring 
with a different geometry (such as a triangle) is used to point towards the two baselines 
simultaneously, all muons will be injected into this ring, but the straight length towards each 
detector will be smaller than in the racetrack case, i.e., 0<SF<2 in general. The scale factor 
is then convenient to parameterize the obtained luminosity relative to the IDS-NF baseline 
setup: SF>1: higher luminosity, SF<1: lower luminosity. Note that, in principle, the SF 
can, for lower E^^, also be increased by a re-optimization of the front-end and generally will 
increase for lower energies due to the reduced decay losses during acceleration. For example, 
a SF=2.8 for a low energy 4 GeV Neutrino Factory has been obtained in Ref. (35] compared 
to SF=2.0. We will not consider this type of effect, since it depends on the accelerator 
complex in a non-trivial fashion. 

For the updated detector simulations, we use the migration matrices mapping the incident to 
the reconstructed neutrino energies for all individual signal and background channels, which 
can be directly implemented into GLoBES. Note that charge mis-identification, (electron) 
flavor mis-identification and neutral current backgrounds are included. For the binning. 



we then follow Ref. 23 , 24 , where the migration matrices for the appearance channels are 
given. For the disappearance channels, we use the same matrices]^ In addition, we increase 
the number of sampling points for high energies to avoid aliasing. This implementation will 
be used throughout the remainder of this paper, unless indicated otherwise. It is denoted 
by the label "new-NF". Note that we also include signal (2.5%) and background (20%) 
normalization errors, uncorrelated among all oscillation channels. 

For the u^- contamination, we use the migration matrix from Ref. |27j for both the z/g — )■ z/^ 

■^That is somewhat on the conservative side, since we require charge identification and better results may 



be obtained with an event sample without charge identification 17 
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Signal 


NC bckg 


CC bckg 


Vt- bckg 


(app) 


7521 


20 


25 


142 


z/^ (app) 


924 


45 


39 


13 


z/^ (disapp) 


4.0 X 10^ 


31 




8154 


z/^ (disapp) 


2.4 X 10^ 


8 




4337 



Table 1: The expected event rates for new-NFr in a 50kt detector at a 4000 km baseline with a muon 
energy of 25 GeV. The chosen oscillation parameters are taken from Eq. Q with 6*13 = 5.6° and (5cp — 0. 



and Vfj, — )■ zz-r channels, since it only depends on characteristics of the r decays. Note that, 
since the binning given in there is different from Refs. 23 ,24 , we had to re-bin this matrix 
carefully. As an important consequence, all events below 2 GeV are collected in the lowest 
bin. We also apply the muon kinematic cuts for the muons from the r decays as for the 
golden channel, following Ref. 



27 



In a more refined approach, one may want to have the 
migration matrices from incident Vr energy to reconstructed energy directly. This setup 
will be denoted as "new-NFr" and it contains everything in new-NF plus the muons from 
r decays. As we will show new-NFr produces practically the same results as new-NF^ 

The input oscillation parameters are taken as follows [2], unless noted otherwise: 



e 



12 



34.4° , e 



13 



5.6° , e 



23 



42.8° 



Am2i 



7.59 X 10"^ eV^ 



2.46 X 10"'^ eV^ 



(7) 



We impose external la errors on Aml^ (4%) and 9i2 (4%) and on Am^i (10%) and 6*23 
(10%) as conservative estimate^ for the current measurement errors |2]. We also include 

Unless noted otherwise, we simulate the normal 



a 2% matter density uncertainty |43||44 
hierarchy. 



2.2 Event rate comparison 

In Fig. [1) we compare the event rates of the latest detector simulation new-NF (thick solid 
curves) with IDS-NF 1.0 (thin solid curves) for the four different oscillation channels as 
given in the plot legend. 

IDS-NF 1.0 (thin curves) did not use any migration matrices and this is reflected in the 
background shape, both neutral current (NC) and charged current (CC), which closely 
follows the signal shape. The signal shape of IDS-NF 1.0 is quite similar to the one of 
new-NF, indicating that migrations are not large for the signal, which is not surprising 
since energy reconstruction works well for the signal events. The background shapes, on 



^This statement is true only for the performance indicators used in this paper, which all focus on the the 
appearance channel, and will most likely not apply to precision measurements of the atmospheric neutrino 



parameters in the disappearance channels as indicated in Ref. 26 . 

^Here we expect that the best measurement of the atmospheric parameters comes from the Neutrino Fac- 
tory in the near future. However, it is useful to add the current information on the atmospheric parameters 
to speed up the marginalization and degeneracy finding. 
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Figure 1: A comparison of the event rate spectra between new-NF 24 (thick curves), including back- 



grounds from i^r 26 27 , and IDS-NF 1.0 (thin curves) for the different oscillation channels as given in the 



plot legend. The chosen oscillation parameters are taken from Eq. ([t]) with ^cp = 0. The muon energy is 
25 GeV and the detector mass is 50kt at a baseline of 4000 km. 
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the other hand, differ substantially between IDS-NF 1.0 and new-NF, since here migrations 
are non-negligible. In particular for the NC background, we observe that for new-NF (thick 
curves) it is quite peaked at low energies. This phenomenon is known as "feed-down" : for a 
given incoming neutrino energy, there will be less energy deposited in the detector in a NC 
event than in a CC event, simply because a neutrino is leaving the detector carrying away a 
sizable fraction of the incoming energy. If a NC event is mis-identified as being a CC event|^ 
then the CC event kinematics will be used for energy reconstruction, which assumes that 
j^vec _ El^^^^^ + -Ehadrons- This results in a systematic downward bias in the reconstructed 
energy for NC background events. This feed-down is the strongest effect of migration and 
thus has potential impact on the energy optimization, since it penalizes neutrino flux at 
high energies, where there is little oscillation but a large increase in fed-down background. 
Also, for muons from r decays there is a strong feed-down for a similar reason: in the decay 
of a r there will be two additional neutrinos which leave the detector. Here, the disruptive 
effect of high energies is even more pronounced, since the Ur CC cross section is a steeply 
increasing function of neutrino energy up to about 30 GeV. 

In summary, the CC backgrounds in new-NF pile-up at lower energies. These low energy 
events are relevant for degeneracy resolution, especially for intermediate values of sin^ 2^13 ~ 
10~^ — 10~^. However, the oscillation peak in vacuum would be at about 10 GeV, and matter 
effects are most important at about 8 GeV, which need to be covered especially for small 
sin^2^i3, where the event rates otherwise rapidly decrease with distance. The backgrounds 
from T decay in new-NFr tend to collect around 8 GeV and may present an immediate 
problem for all values of sin^ 26'i3. Therefore, it is not quite clear that high muon energies 
are preferred everywhere in the parameter space, and one may suspect that the baseline- 
muon energy optimization may be a complicated function of the detector response. 



2.3 Performance and impact of Vt- contamination 

Neutral current backgrounds do not carry any information about flavor conversions of active 
neutrinos and therefore are detrimental to oscillation searches. The muons from r decays, 
on the other hand, do arise from oscillation and they are a sign of appearance of a new 
flavor, r, in a beam otherwise devoid of this flavor. The background arising from as 
defined in Eq. ([T]) (appearance channels) and Eq. ^ (disappearance) are shown as gray 
(brown) solid curves in Fig. [l] In all channels, they are the largest source of background. 
It is, however, not clear from the beginning whether this is a benefit or a curse, since 
this oscillating background carries information on the oscillation parameters. In particular, 
the low energy parts, which actually stem from much higher incident neutrino energies, 
may carry complementary information to the high energy signal; since the resulting energy 
distribution is different they may be separated on a statistical basis. For example, the z/^ 
appearance probability is given, expanded to second order in sin26'i3 and the hierarchy 



parameter a = Aml^/ Aml-^ ~ 0.03, as |13 45[46| : 
- sin^ 26*13 sin^ 6*23- 



sin2[(l-i)A3i] 



^Otherwise, it would not be a background event. 
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Figure 2: A comparison of the discovery reach of CPV, MH, and 0i3 at the 3cr CL among different 
detector simulations. The label "IDS-NF" refers to the detector in the IDS-NF baseline setup 1.0 The 



simulation with the migration matrices from 23 is indicated by the label "1004.0358". The label "new- 



NF" refers to most up-to-date detector simulation in Ref. p4]. The Vr contaminations in the appearance 



and disappearance channels are, in addition, included in "new-NFr" 26 27 . Here a combination of two 



baselines 4 000 km and 7 500 km with two 50 kt MIND detectors is assumed. 



± a sin 29i3 sin 5cp sin 29i2 sin 26*23 sin ( A31 ) 
+ a sin 26*13 cos 6cp sin 26^12 sin 26*23 cos(A3i' 
+ c? cos^ 6*23 sin^ 26*12 



sin(iA3i) sin[(l - i)A3i] 
sin(iA3i) sin[(l - i)A3i] 



A 



[\-A) 



.2 2n ..„2ozi sin2(iA3i) 



with A31 = Aml^L/{AE) and A = ±2 y/2 E Gp nJAml-^. The signs in the second term 
and A are positive for neutrinos and negative for anti- neutrinos. For Pgr, the channel which 
controls the background in Eq. Q, flip the sign of the second and third terms and replace 
in the first and fourth terms sin^ ^23 ^ cos^ 6*23. For maximal atmospheric mixing, only the 
signs of the second and third terms change. Now consider, for instance, the magic baseline 
L ~ 7500 km where, by definition, sin(AA3i) ~ 
survives, which is the same for the signal and for the 
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In this case, only the first term 



background, which means that it adds 
to the sin^ 2^13 and MH sensitivity. For the short baseline used for the CPV measurement, 
the sign of the second and third terms are different between Pg/x cind Per, which means 
that the effects of 6cp are, naively, reduced by the z/^ background. However, note that 
the background is reconstructed at lower energies, which means that one can, in principle, 
distinguish the two channels. It is therefore, without numerical simulation, not obvious if 
the Ur contaminations improve or deteriorate the sensitivities. 

The physics performance arising from the different detector simulations for the CPV, MH, 
and 6^13 discovery reaches are show in Fig. [2] Here "IDS-NF 1.0" refers to the detector 
performance of the IDS-NF baseline setup 1.0 [5]. The results in the figure demonstrate that 
the performance based on the detector simulation presented in Ref. [23] (thick solid curves) 
is worse. The main reason, we were able to identify, is significantly higher backgrounds from 
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Figure 3: A comparison of the discovery reach of CPV, MH, and 6*13 at the 5a CL among different 
experimental setups: "50kt+50kt" refers to a combination of 50 kt MIND at 4 000 km and 50 kt MIND at 
7 500 km (SF=1), "lOOkt only" to a 100 kt MIND at 4000 km (SF=2), "100kt+50kt" to the combination of 
a 100 kt MIND at 4 000 km and 50 kt MIND at 7 500 km (SF=1). All these setups use the most up-to-date 
detector simulation new-NF |24 . "IDS-NF 1.0" refers to IDS-NF baseline setup 1.0, i.e., the combination 
of 50 kt MIND at 4 000 km and 50 kt MIND at 7 500 km (SF=1), using no migration matrices [s] (note 
number IDS-NF-002), to be compared to the dotted curves. 



charge mis-identification than in the IDS-NF 1.0. The most up-to-date detector simulation 
is presented in Ref. [24] (dotted curves) and this setup is labeled new-NF, for which the 
performance is slightly better than for the IDS-NF 1.0. In this case, the signal efficiencies and 
threshold are improved compared to IDS-NF 1.0, while the background level is maintained. 
One of the main differences with respect to Ref. 23 is the inclusion of quasi-elastic events 
which improves the signal efficiency at low energies. The effect of the migration of the 
backgrounds does not have a large impact on the discovery reaches. This may not be true 
for precision studies of the atmospheric oscillation parameters, however, a detailed answer to 
this question is beyond the scope of the current paper. If, in addition, the contributions from 
the Ut- are included, new-NFr (dashed curves), there is hardly any effect on the performances. 
Note, that the relative impact of r decays does depend on the underlying detector parameters 
and this illustrates that it is difficult to predict the effect of the z/^ without numerical 
simulation. In any case, the absence of a significant difference in performances between 
new-NF and new-NFr is in agreement with the results presented in Ref. 27 and therefore, 
we will not further consider r decays and the resulting backgrounds. 

Other questions to be addressed in the context of the updated detector simulation are the 
quantitative comparison between one and two baselines, and the impact of a larger detector 
at the shorter baseline. We discuss these in Fig. [3} where several versions of the updated 
detector are compared with the IDS-NF 1.0. Note that the scale factor (SF) has been 
adjusted for the assumed racetrack storage rings to correct for the larger number of use- 
ful muon decays during the single baseline operation. In addition, note that this figure is 
shown at the 5a CL, compared to the previous, to make the impact of degeneracies clearer. 
From the comparison of the IDS-NF 1.0 and the corresponding 50kt-|-50kt curves using 
new-NF confirm the earlier result, note that there is not much difference in performance. 
A possible alternative setup is to operate a single 100 kt detector at the 4000 km baseline. 
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this configuration is labeled "lOOkt only" . This configuration actually exhibits better per- 
formances for CPV and the ^13 discovery because of the factor of two higher luminosity 
using the racetrack-shaped storage rings. In this case, the complementary information at 
the 7500 km baseline is replaced by high statistics at the short baseline. Note, however, 
that the MH discovery reach is significantly worse, and that degeneracies affect the shape 
of the CPV curve. The setup "lOOkt-FSOkt" , where there is a 100 kt detector at 4000 km 
and one 50kt at 7500 km, can easily resolve the degeneracies at about sin^ 2^13 ~ 10~^ in 
the CPV discovery reach, while the MH and ^13 discovery reaches are comparable. In this 
case, SF=1, which means that this setup in fact only has 75% of the exposure of the "lOOkt 



only" version. Therefore, the two baselines are synergistic in the sense of Ref. 47 , i.e., for 



the same exposure, the baseline combination clearly performs better. However, if one sticks 
to the racetrack geometry of the storage rings, the one baseline operation may be more 
efficient. For a triangular shaped ring, which we will discuss later, this argument changes, 
because the second baseline is available anyway. The question of the necessity of the magic 
baseline remains open. Especially in the context of new physics and surprises, such as a 
lower than expected machine luminosity, it provides a robust alternative. 

In the following, now that we have quantified the impact on the performance, we will only 



consider the setup with the updated migration matrices from Ref. 24 , i.e., new-NF. Wher- 



ever we refer to "IDS-NF" , we will actually mean the IDS-NF parameters {E^ = 25 GeV, 
4000 km + 7500 km), while the detector simulation is new-NF. We will not consider the 
Ur contribution anymore, partially because it has been shown not to have a significant im- 
pact on the discussed performance indicators, partially because it will need to be quantified 
within the same detector simulation as the signal and other backgrounds in the future. 



3 Optimization of a green-field setup, low versus high 
energy Neutrino Factory? 

Here we study the optimization of a green-field setup, which means that no particular 
accelerator and detector sites are chosen and that the baselines and muon energy are not 



constrained. The optimization is performed using the migration matrices from Ref. 24 
Now that the detection threshold has improved, we are especially interested if the new 
MIND detector can interpolate between low and high energy Neutrino Factory. 

First of all, consider that sin^ 2^13 is not found before the Neutrino Factory operation. 
Assume that, in this case, one wanted to optimize for the reach in sin^2^i3, i.e., CPV, 
MH, and 613 should be discovered for as small as possible true values of sin^2^i3. For 
the sake of simplicity, we choose maximal CP violation 6cp = 7r/2 for the true 5cpQ We 
show in Fig. the discovery reach in sin^ 26'i3 for maximal CP violation, MH, and ^13 as 



a function of baseline and E^. The contours show the reach in (true) sin 2^13 for which 
the different quantities will be discovered at the 3a CL. This figure is to be compared 



^Other, more technical versions, are choosing the "typical value of (5cp" (the median of the distribution 
in 6cp), corresponding to a fraction of 6cp of 50%, or a different certain fraction of 6cp- At least for CPV, 
our choice corresponds to the most optimistic case. 
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Figure 4: Discovery reach in sin^ 2^13 for maximal CP violation, JVIH, and ^13 as a function of baseline 
and E^. The contours show for how small (true) sin^ 26*13 the different quantities will be discovered at the 
3cr CL, where maximal CP violation ^cp = 7''/2 is chosen as a true value in all cases. The best reaches for 
basehne and are marked by dots: (4519,16.25), (5805,22.57) and (4800,22.50) followed by their optimal 
sensitivity of sin^ 26*13 at 10"'*•^ lO""-^ and lO-'* ^ Here SF=1 is used with one 50 kt detector. 



to Figs. 5 and 6 of Ref. 17 for the respective ^cp and an older version of the detector 
simulation. Here the qualitative features are clearly recovered: The CPV discovery requires 
a 2 500 km to 5 000 km baseline and above about 12 GeV. Note that degeneracies are 
typically unproblematic for this choice of 5cp, whereas for 5cp = 37r/2, a second baseline 
may be required. In addition, note that relatively low are allowed because of the low 
detection threshold. For the MH discovery, baselines longer than 4 000 km and i?^ larger 
than about 10-12 GeV are needed, since the MSW resonance energy of about 8 GeV is to 
be covered. Here even longer baselines are preferred for different values of 5cp- For the 
discovery, we find an extremely wide baseline and energy range, giving the least constraints. 
However, note again that this result depends on the choice of 5cp- In summary, the result 
of this optimization, qualitatively, points towards one baseline between 2 500 and 5 000 km 
for the CPV measurement and one very long baseline for the MH measurement, such as 



the magic baseline at 7500 km useful for degeneracy resolution (see Ref. 17 for a more 
detailed discussion). Because of the optimized detector, lower E^ of down to 12 GeV may 
be possible. Below, we will discuss how this result changes for specific true values of sin^ 26i^ 
if all values of 5cp are considered. 

From a different perspective, consider that the value of sin^ 26'i3 is known, either from an 
earlier stage experiment or an earlier stage of the Neutrino Factory. In this case, as we 
have seen in the previous section, the MH discovery is typically not a problem (at least in 
combination with a longer baseline if sin^ 26*13 is small), and the most interesting question 
is the optimization of the fraction of 5cp for which CPV can be discovered. We first show 
in Fig. [5] the fraction of 5cp for which CPV will be discovered (Sa CL) as a function of 
L and E^ for the single baseline Neutrino Factory. The different panels correspond to 
different true values of sin^ 26*13, as given there. From this figure, it is obvious that the 
optimization strongly depends on the value of sin^ 2^13 chosen. For large sin^ 26*13 — 10~^, 
shorter baselines and lower energies are preferred. Even E^j, as low as 5 GeV at the FNAL- 
Homestake baseline of about 1 300 km is not far from optimal, which means that the MIND 
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Figure 5: Fraction of Sqp for which CPV wiU be discovered (3ct CL) as a function of L and Ef^ for 
the single baseUne Neutrino Factory. The different panels correspond to different true values of sin'^ 2(?i3, 
as given there. Here SF=1 is used with a 50 kt detector. The optimal performance is marked by a dot: 
(2200,10.00), (2288,13.62), (3390,20.00) and (4345,22.08) with regard to their best reaches of the fraction 
of ^cp at: 0.77, 0.84, 0.67 and 0.42. 
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detector approaches the TASD performance of the low energy Neutrino Factory. Very 
interestingly, compared to earlier analyses without background migration, too high E'^ are 
in fact disfavored in the large sin^ 2^13 case. Note that neither for the considered detector 



nor for the TASD, we find strong evidence supporting the "bi-magic baseline" in Ref. 48 
see discussion in Appendix [cj For the other extreme, sin^ 29i3 — 10~^, baselines between 
4 000 and 5 000 km are preferred with ii^^ ~ 20 — 25 GeV, which corresponds more to 
the high energy Neutrino Factory, such as the IDS-NF baseline. Including the other two 
panels, the optimal region within each panel moves from the lower left on the plots to the 
upper right as the value of sin^ 2^13 decreases. This means that, depending on the choice 
of sin^2^i3, the optimization results in the low energy Neutrino Factory, the high energy 
Neutrino Factory, or an intermediate scenario, and that the low and high energy Neutrino 
Factories are just two versions of the same experiment in different optimization regions. 
Of course, this discussion is somewhat hypothetical from the practical point of view, since 
either the next generation(s) of experiments will find sin^ 2^13 or not. If they find sin^ 2^13, 
the optimal parameters of the Neutrino Factory can be clearly predicted as a function of the 
detector response. The FNAL-Homestake low energy Neutrino Factory is one such possible 
setup for large enough sin^ 2^13 for the MIND detector. If they do not find sin^ 26'i3, one 
may want to go for the IDS-NF high energy setup, which, in a way, represents the most 
aggressive but also inclusive option: This version of the Neutrino Factory is optimized for 
the worst case scenario. 

Apart from the single baseline, we show in Fig. [6] the combination with another fixed base- 
line L2 = 7500 km (in fact, there is typically very little dependence on the exact choice 
of the second baseline [iS]). Note that the muon energy is the same for both baselines. 
Comparing Fig. |6] with Fig. |5| we find that the optimization of the short baseline hardly 
changes for very small and very large sin^2^i3, whereas the possible baseline windows for 
intermediate sin^ 2^13 (upper right and lower left panels) become somewhat broader. The 
energy optimization remains almost unaffected. As far as the absolute performance is con- 
cerned, especially for sin^ 26*13 = 10~^ and sin^ 26*13 = lO""^, the fraction of ^cp increases 
because of the degeneracy resolution potential of the second baseline (which is not sensitive 
to Sep itself by choosing exactly the magic baseline). For large values of sin^ 26*13, the sec- 
ond baseline is not required. This again reflects the correspondence to low and high energy 
Neutrino Factory: the low energy version is typically proposed with one baseline, the high 
energy version with two baselines. 



4 Earth geometry, and triangular shaped storage ring? 

In this section, we discuss the geometry aspects of specific sites for the high energy Neutrino 
Factory. The relevant questions for us are: 

1. Can we find possible baseline combinations for the high energy Neutrino Factory for 
the large accelerator laboratories on different continents? 

2. Would it be possible to use a single, triangular-shaped storage ring pointing towards 
both detector locations at the same time? 
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Figure 6: Fraction of Sep for which CPV will be discovered (3cr CL) as a function of Li and Ef^ for the 
two-baseline Neutrino Factory, where L2 — 7 500 km fixed. The different panels correspond to different true 
values of sin^ 2^13, as given there. Here SF=1 is used and E'^ is assumed to be equal for both baselines. 
Two 50 kt detectors are used in the simulations. The optimal performance is marked by a dot: (2300,10.00), 
(2580,15.00), (3250,20.00) and (4297,22.05) followed by their best discovery reach of the fraction of 6cp at 
0.77, 0.84, 0.81 and 0.51. 
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CERN 


FNAL 


J-PARC 


RAL 




(46.24,6.05) 


(41.85,-88.28) 


(36.47,140.57) 


(51.57,-1.32) 


Asia: 










CJPL (28.15,101.71) 


7660 


10420 


3690 


7840 


Kamioka (36.14,137.24) 


8770 


9160 


300 


8640 


YangYang (37.77,128.89) 


8350 


9300 


1050 


8270 


INO (9.92,78.12) 


7360 


11410 


6570 


7820 


Europe: 










LNGS (42.37,13.44) 


730 


7350 


8840 


1510 


Pyhasalmi (63.68,25.98) 


2290 


6630 


7090 


2080 


Slanic (45.27,25.95) 


1540 


7780 


8150 


2110 


Boulby (54.56,-0.81) 


1050 


5980 


8480 


340 


Canfranc (42.76,-0.51) 


650 


6550 


9280 


980 


Frejus (45.20,6.67) 


130 


6830 


8900 


920 


SUNLAB (51.22,16.16) 


930 


6980 


8190 


1210 


Umbria (42.98,12.64) 


640 


7280 


8830 


1420 


Gran Canaria (28.39,-16.59) 


2780 


6240 


10570 


2850 


North America: 










Soudan (47.82,-92.24) 


6590 


730 


8500 


5900 


WIPP (32.37,-104.23) 


8160 


1760 


8900 


7540 


Homestake (44.35,-103.77) 


7360 


1290 


8250 


6690 


SNOLAB (46.47,-81.19) 


6090 


760 


8950 


5400 


Henderson (39.77,-105.86) 


7750 


1500 


8410 


7110 


Icicle Creek (47.56,-120.78) 


7810 


2610 


7240 


7160 


San Jacinto (33.86,-116.56) 


8600 


2610 


8170 


8000 


Kimballton (37.37,-80.67) 


6580 


820 


9560 


5950 



Table 2: Here we show the basehnes between the considered accelerator facihties (columns) and under- 
ground laboratories (rows) in kilometers. The latitude and longitude of each site is given in the brackets 
in degrees. The baselines are calculated using Mathematica with the International Terrestrial Reference 



Frame 2000, rounded to 10 km. All coordinates are consistently extracted from Google Maps 49 



We will quantify in the next section how specific baseline combinations translate into per- 
formance and optimization compared to the IDS-NF baseline parameters. 

In order to address these purely geometric questions, we consider CERN, FNAL, J-PARC, 
and RAL as potential host laboratories for the Neutrino Factory]^ For the potential detector 
sites, we adopt the conservative point of view that significant rock overburden is needed. 
This assumptions and the anticipated timescale of the Neutrino Factory limits the choice 
of potential detector sites to currently investigated, or at least discussed, deep underground 
laboratories. We list the potential accelerator facilities and underground laboratories to- 
gether with their locations and baselines between them in Table |2j see Appendix |A] for more 
details on the individual locations. The locations of laboratories and detector sites on the 

®Note that J-PARC is not very far away from KEK, for which the discussion would hardly change. 
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Figure 7: A world map including all potential sites of accelerators with large symbols and underground 
labs with stars. 

Earth's surface can be found in Fig. [7j 

For CPV, the IDS-NF basehne has been 3000 km to 5000 km, based on the analysis in 
Ref. [It] . This conclusion was obtained from the optimization of the ^13 reach, similar to 
Fig. |4|, left panel. However, from Fig. [s] and Fig. |6] (see also Fig. 8 in Ref. [37]), those shorter 
baselines are preferable if sin^ 2^13 turns out to be somewhat larger. Therefore, we allow for 
Li G (1500, 5000) km for the high energy Neutrino Factory. For degeneracy resolution and 
the mass hierarchy measurement, L2 should be close to the "magic baseline" [iTj (Fig. Elis 
for one specific true value value of ^cp); which can be see, for instance, in Fig. 5 of Ref. p!9] . 
This location does not need to be exact. However, the baseline should not be to short, in 
order to allow matter effects to pile up and to suppress the CP violating terms, and not too 
long if too steep active storage rings legs should be avoided. We choose L2 G (7000, 8000) 
km clS cl reasonable range, see Ref. jl8j. 

As we can read off from Table [2| there is a very limited number of the short-baseline Li 
detector sites: 

CERN Li. Pyhasalmi (Finland), Slanic-Prahova (Romania), Gran Canaria (Spain). 
FNAL Li'. WIPP, Henderson, Icicle Creek, San Jacinto. 
J-PARC Li: CJPL. 

RAL Li. LNGS (Italy), Pyhasalmi (Finland), Slanic-Prahova (Romania), Gran Canaria 
(Spain). 
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Fi gure 8: Baseline combinations from Table [s] shown in a two baseline plot. The different colors/symbols 
represent the different laboratories. 



We have not found any baseline between 4 000 and 5 000 km. Obviously, we have plenty of 
options on the second baseline in comparison with the number of the first baseline. It may 
be noteworthy that CERN-INO and CERN-Homestake are exactly the same distances. 

Now with these baseline windows for the short and long baselines, we can, for each lab- 
oratory, choose all possible combinations from Table [2j We show these in Fig. |8} with 
different shapes and colors for each laboratory, and we list them in Table [3j Note that 
several qualitatively different baseline combinations in Fig. [s] are marked (with the numbers 
from Table [3]), which we will discuss in the next section. In addition to the setups with 
the above criteria, we have listed one option with FNAL-Homestake as first baseline (#51). 
As we will demonstrate later, this baseline may be too short for the high energy Neutrino 
Factory. 

Depending on the two baseline combination, it may be possible to use a triangular shaped 



storage ring instead of two racetracks. Here we follow the discussion in Ref. 10 , which the 
IDS-NF baseline setup with two storage rings is based on. The two racetrack-shaped storage 
rings are assumed to have a circumference of 1609 m. The active straights are about 600 m 
long, and, in each storage ring, /z"*" and fi~ circulate in different directions. For a triangular 
shaped ring, probably two beam lines in the same tunnel are required to store /i+ and /i~ 
simultaneously. We assume that the circumference of the triangular ring, representative 
for the tunneling cost, is the same as for one racetrack, and we assume a (conservative) 
curvature radius Rc of about 78 m for the curved sections. For the sake of simplicity, we 
only consider isosceles triangles with the same useful number of muon decays for the two 
far detectors. In the racetrack design of the IDS-NF baseline setup, 2.5 ■ 10^° useful muons 
per year and polarity decay in each straight of each storage ring. For a triangular ring, 
all muons of one polarity can be injected into the same ring, leading to 5.0 ■ lO^'' useful 
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Table 3: Considered two-baseline combinations (sec main text). The five right columns give the pa- 
rameters of an isosceles triangle as storage ring: straight length (meters), dead length (meters), apex angle 
(degrees), scale factor SF (compared to two racetracks), angle to vertical (degrees). 
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Figure 9: Examples for three isosceles triangular-shaped rings (maximum efficiency, racetrack-like 
efficiency, and minimum efficiency, respectively). The numbers refer to Table [sj 

muon decays per year and polarity over a straight length of 600 m, which corresponds to 
SF=2. Of course, due to the fixed circumference, the straights will be shorter than 600 m 
for the triangle, i.e., SF<2. Note that SF=1 corresponds to the same luminosity of the 
two designs. Obviously, if SF>1, the triangle is more efficient, possibly with a factor of two 
lower tunneling cost (because only one tunnel is needed). If SF<1, the loss of efficiency 
could be compensated by a slightly larger storage ring. For SF<0.5, the two racetracks 
are definitively the better option. Note that synchtrotron losses in the curved sections of 
the storage rings and other beam losses are taken into account for the racetrack-shaped 
geometry, yielding the anticipated number of useful muon decays for the IDS-NF baseline. 
We do not expect that these losses change drastically for a triangular-shaped ring. For 
example, for the curved sections, we have assumed the same curvature radius, and the three 
curves sections add up to a circle just as for the two curves sections of a racetrack-shaped 
ring. Since the circumference of the triangular and racetrack shaped rings is the same, the 
muons are exposed to the same curved section distance over their lifetime. 

As the first observation, a triangular shaped ring can always be built if the circumference of 
the ring is larger than 271 Rc {Re- curvature radius curved section), which we have satisfied. 
This can be easily seen by the fact that in the smallest (extreme) case, the triangle with 
curved sections will collapse into a circle (with zero straight lengths). For larger triangles, 
the efficiency of the active legs (SF) may still be extremely small. However, we list in Table|3] 
the triangular geometry in terms of straight lengths, dead section length, apex angle, SF, and 
V-angle (angle between triangle plane and vertical), and it turns out that 0.95 < SF < 1.38. 
In the optimal case (#30), the V-angle is 90°. This means that a triangle could be built for 
all of the considered options. We show three examples for maximum efficiency, racetrack- 
like efficiency, and minimum efficiency, respectively, in Fig. |9j where also the numbers from 
Table [3] are given. In the extreme cases, the triangle resembles a racetrack with either a 
very short or very long dead section. In the worst case, if the two detector locations are 
quite aligned, less then 50% of the useful muon decays over the whole ring can be used. 
However, the factor of two higher muon injection rate compensates for that. 



19 



In this discussion we have ignored how deep the tunnels would be and that two racetracks 
have other advantages. For instance, if one racetrack or one detector needs maintenance, 
all muons can be injected into the other storage ring without loss of performance integrated 
over the whole operation time. However, this discussion is interesting from a different 



perspective: Earlier in Sec. |2.3[ we have shown that a single baseline operation may be 
more beneficial in parts of the parameter space, where one of the reasons is a factor of two 
gain in exposure compared to the operation of two racetracks. However, if a triangular 
ring is built, the argument changes. In Fig. [3} the 100kt+50kt option has an exposure of 
150 kt*l (SF)=150 kt and the lOOkt option an exposure of 100 kt*2 (SF)=200 kt . For 
the triangle, one in the most optimistic case for the 100kt+50kt option has an exposure of 
150 kt*1.38 (SF)~200 kt, which is the same as for the one baseline case - at a much better 
sensitivity, and with the same storage ring circumference. From a different point of view, 
one has the performance depicted by 100kt+50kt in Fig. [3] in that case already with two 
72 kt and 36 kt detectors. Thus, for small sin^ 2^13 and proper detector sites, the triangle 
may finally be the better choice. Note that in the following, unless noted otherwise, we do 
not use the SF from Table [3| but use SF=1 instead (two racetracks). 

In summary, we have demonstrated that reasonable pairs of detector locations can be found 
for the considered accelerator laboratories. We have stated that one could always use 
a triangular-shaped storage ring with a similar efficiency as two racetracks from purely 
geometrical arguments, and that the efficiency varies at about 40% among the different 
options. 



5 Site-specific performance and energy optimization 

Here we discuss the performance of site-specific setups, as well as the optimization of E'^ for 
specific sites. Because of the large number of options considered, we only show examples in 
this section, whereas the plots for all discussed sites can be found in Appendix [B] 

Let us first of all quantify the performance in comparison to the IDS-NF baseline combi- 



nation 4000 km+7500 km at SF=1. Therefore, we show in Fig. 10 the discovery reach for 
CPV, MH, and 6^13 discovery {3a) for a number of qualitatively different selected basehne 
combinations for different accelerator laboratories (in rows). For each laboratory, we have 
chosen an example roughly representing the best case and an example close to the worst case 
for the chosen = 25 GeV, as well as we show the IDS reference values. Here two racetrack- 
shaped storage rings (SF=1) are assumed. In all accelerator cases for CPV, options can 
be found which perform better than the IDS combination if sin^ 2^13 > 10~^, because large 
values of sin^ 26'i3 prefer shorter CPV baselines, as discussed earlier. In these cases, even 
a single baseline option with a lower could be preferable. For 10^^ < sin^ 26^13 < 10~^, 
options close to the IDS performance can be easily identified. For sin^ 2^13 ^ 10~^, the IDS 
combination can roughly be matched, but the sensitivity cannot be exceeded, at least not 
with the racetrack-shaped storage rings. The reason is that we do not use any baselines close 
to, or exceeding 4 000 km. Because of the absence of potential detector sites, one may want 
to study either alternative locations, or the possibility to use MIND close to the surface. In 
this case, the long baseline may actually help for background suppression, since neutrinos 
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sin^26'i3 sin^26'i3 sin^2(9i3 



Figure 10: The discovery reach of CPV, MH, and ^13 (3tT, m columns) for selected baseline combinations 
for different accelerator laboratories (in rows). In all panels, the curves for the IDS reference combination 
4 000 km+7 500 km with new-NF is shown for reference. Here i?^ is fixed to 25 GeV with two 50 kt detectors, 
and SF=1 is used in all cases (two racetrack-shaped storage rings). 
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from directions close to the beam have to travel through a significant amount of rock then. 
In the worst case scenarios, significant sensitivity losses may have to be taken into account, 
especially if not long enough CPV baselines are used. The MH discovery, on the other 
hand, is driven by the long baseline, but also benefits from a longer short baseline. For the 
optimal options, the IDS sensitivity can be matched, although it may be a bit different as 
a function of Sqp- Similar results are obtained for the 613 discovery reach. 

In Appendix [B| we show the performances of the discussed combinations (#1 to #50) from 
Tablejs} In each case, we compare the performances with racetracks (SF=1, dotted curves), 
triangular-shaped geometry (SF from Table |3| dashed curves), and IDS combination (SF=1, 
solid curves). Here we just mention some of the most interesting options from these figures, 
especially those with excellent sensitivity for sin^ 29i3 ^ 10~^ which can be further improved 
by a triangular shaped ring. Here CERN or RAL to Gran Canaria and to CJPL or INO 
are interesting options with a significant sensitivity gain and good absolute performance. 
In addition, J-PARC to CJPL and Icicle Creek is, in fact, the only option we find which 
can exceed the IDS reference performance for small sin^2^i3. For options with shorter 
baselines, the performance also improves significantly in many cases by using a triangular 
geometry, but that cannot compensate for the baseline choice. In no case, the performance 
is significantly worse using a triangle. 

Another interesting question is the optimization of for specific two baseline setups. 
Remember that the two baseline Neutrino Factory is mostly relevant for small sin^26'i3. 
Here we choose the combinations from RAL as example, see Fig. 11, the other laboratories 
are shown in Appendix |Bj Basically, we can identify three different sets of curves in that 
figure, which correspond to the three different CPV baselines in Table [s] (a) Li = I 510 km, 

(b) Li ~ 2100 km (two different ones), (c) Li = 2 848 km. Depending on the value of 
sin^ 26*13 and i?^, one of these three sets performs best: below about 7-8 GeV (depending 
on sin^ 26^13), (a) is best, between about 8 and 12 GeV (b) is best, and above 10-14 GeV, 

(c) is best. This results more or less reproduces the green-field optimization. Note that for 
sin^ 2^13 < 10~^, the case for which the two baseline Neutrino Factory is the relevant choice, 
the long CPV baseline options are better in terms of absolute performance, provided that 

is high enough. In addition, note that the IDS reference prefers > 20 GeV in all 
cases, where the discovery reaches saturate. To summarize, the optimal muon energy does 
not only depend on sin^ 2^13, but also on the specific two baseline combination. However, in 
many cases, the performance saturates at about 12-15 GeV (see, e.g., sin^ 26*13 = 10"'^), and 
in some cases may even decrease for too high E^. The IDS-NF baseline choice Efj_ = 25 GeV 
can be understood as an aggressive option from the current point of view. However, note 
that for any given baseline combination, the optimization of E^ can be easily performed. 
From the machine point of view, it should be easy to "down-grade" the setup then. 



6 Summary and conclusions 

In this study, we have revisited the optimization of the Neutrino Factory based on the 
most up-to-date analysis of the MIND detector using migration matrices. We have also 
considered possible backgrounds from taus, which come from charged current interactions 
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CPV@RAL 




[GeV] [GeV] 

Figure 1 1 : CPV discovery reach at RAL as a function of E'^ for all RAL-options in Table [s] at the 
3ct CL for different values of (true) sin^ 2^13 (as given in the panels). Here we assume SF=1.0 and two 
50 kt detectors. We group the different baseline combinations according to the shorter baseline Li, as 
shown in the legend. The bands basically reflect the variation of the second baseline. The IDS-NF baseline 
combination is shown by the black curves with Li = 4000 km and L2 = 7500 km. 




sin 2013 sin 2013 sin 20i3 



Figure 12: A comparison of the discovery reaches of CPV, MH, and ^13 at the Str CL between a low 
energy single baseline neutrino factory (LENF) with E,, ^ 10 GeV and L = 2000km (c/., Fig.js) upper left 
panel), and a high energy two baseline neutrino factory (HENF) with i?^ = 25 GeV, L\ = 4000km, and 
L2 — 7500km (c/.. Fig. [6] lower right panel). 
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in the detector, and which practically cannot be distinguished from muons. We have found 
that the resulting backgrounds do not have a visible impact on the CPV, MH, and 613 
discovery reaches. A more refined discussion will require, however, a consistent treatment 
of all migration matrices. 

Although the optimization of the Neutrino Factory does generically not change with the 
new detector simulation, there are a number of interesting observations. The lower thresh- 
old and higher efficiencies compared to earlier simulations imply that the MIND detector 
characteristics are getting more similar to the characteristics of the detectors proposed for 
the low energy Neutrino Factory {e.g., a magnetized TASD). We could demonstrate that 
we recover the L-ii^^-optimization of the low energy Neutrino Factory for large sin^ 2^13: In 
this case, a single baseline Neutrino Factory with i?^ as low as 5 GeV and a baseline as 
short as FNAL-Homestake (about 1 300 km) might be sufficient. For small sin^ 2^13 < 10~^, 
however, we find that a two baseline Neutrino Factory with one baseline between about 
2 500 km and 5 000 km, the other one at about the magic baseline 7 500 km is ideal. This 



recovers the results from Ref. 17 . We summarize this in Fig. 12, where we compare the 
performance of the optimal single baseline low energy Neutrino Factory with the optimal 
two baseline high energy Neutrino Factory for the same MIND detector. One can clearly 
see that for sin^ 2^13 > 10~^ the low energy version can perform all of the required measure- 
ments, whereas for smaller values the high energy Neutrino Factory is clearly better. This 
means that this different optimization would be sufficient to compensate for the relative 
deterioration of performance at large sin^ 2^13 observed in the traditional high-energy Neu- 
trino Factory. If sin^ 2^13 was known, the shorter (CPV) baseline could even be optimized: 
The larger sin^ 26'i3 was, the shorter CPV baselines would be preferred in the mentioned 
baseline window. The next generation of experiments will tell us if sin^ 26'i3 > 10~^ or 



smaller, see, e.g., Ref. 50 , therefore, we can optimize for the large sin^ 2^13 > 0.01 case. 



Note that a more refined optimization depending on the size of sin 2^13 may be possible 



for a staged Neutrino Factory approach, as it is illustrated in Ref. 37 . 

Apart from the optimization of the green-field Neutrino Factory, we have performed a site- 
specific analysis for the high energy Neutrino Factory, assuming that it requires two base- 
lines. We have considered four different accelerator laboratories on three different continents 
(CERN, FNAL, J-PARC, RAL) and a number of potential detector locations in suggested 
underground laboratories. We have found that in all cases plausible baseline combinations 
can be found. However, for small sin^2^i3, where a baseline between 2 500 and 5 000 km 
is preferred for CPV, we only found one possible baseline: J-PARC to CJPL (China). 
Therefore, we propose that possible underground sites for this baseline window should be 
investigated. In addition, we propose to study the MIND performance on the surface, since 
surface operation would greatly facilitate site and basehne selection. 

We have also investigated the possibility to use a triangular-shaped muon storage ring 
compared to two racetracks, where the efficiency is a function of the Earth geometry and 
the chosen source and detector locations. We have first of all shown that solely based on 
geometry a triangular ring could be used in either case, without significant loss of luminosity. 
Then we have identified a number of baseline combinations with reasonable baseline lengths 
for which the triangle would be especially interesting: CERN or RAL to Gran Canaria and 
to CJPL or INO. In addition, J-PARC to CJPL and Icicle Creek is, in fact, the only option 
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we found which can exceed the IDS reference performance for small sin^ 2^13. We have also 
pointed out that using a triangular-shaped ring, the decision between one and two baselines 
does not emerge, and that, from the physics point of view, the two baseline combination is 
more efficient for the same exposure. 

As far as the optimization is concerned, the feature of the new detector simulation 
that the backgrounds are typically reconstructed at lower energies and that the threshold is 
lower leads to new insights. Especially for the high energy Neutrino Factory in the context 
of specific baseline combinations, the CPV performance in some cases saturates at lower 
muon energies than 25 GeV. Although the current IDS-NF baseline setup with = 25 GeV 
is still optimal for the optimal baseline combination, and high E^ typically do not harm for 
small sin^2^i3, smaller energies may be preferred for specific sites. 

We conclude that the low energy and high energy Neutrino Factory should not be regarded 
as separate options. Let us emphasize that the optimization, for instance, of E^^ is a function 
of sin^ 26*13, the detector response, and the specific sites chosen for detector and accelerators. 
Therefore, the IDS-NF baseline with E^ = 25 GeV should be understood as most conser- 
vative choice, which can be downgraded in specific scenarios/for specific detectors. From 
the machine point of view, we recommend to choose splitting points between the different 
accelerator components at about 5 and 12 GeV, which will allow for E^ = 5, 12, or 25 GeV. 
The final choice has to be made based on the knowledge on sin^ 26*13 at the time of decision, 
the choice of the detector, and the specific site. 
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A Locations of accelerator facilities and underground 
laboratories 



We use Google Maps 49 to find the exact locations (latitudes and longitudes) of considered 
accelerator facilities and underground laboratories. In the following, we provide the details 
of the locations based on which the latitudes and longitudes have been obtained. 
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A.l Accelerator facilities 



• CERN 

Latitude : 46.24° N & Longitude : 6.05° E 

Route de Meyrin 385, 1217 Geneve, Schweiz, Switzerland. 

• FNAL 

Latitude : 41.85° N & Longitude : 88.28° W 

Center for Particle Astrophysics, Fermi National Accelerator Laboratory, P.O. Box 
500, Kirk Road and Pine Street, Batavia, Illinois 60510-0500 USA. 

• J-PARC 

Latitude : 36.47° N & Longitude : 140.57° E 

Tokai Village, Naka District, Ibaraki Prefecture, Japan. 

• RAL 

Latitude : 51.57° N & Longitude : 1.32° W 

Rutherford Appleton Laboratory, Harwell Science & Innovation Campus, Didcot OXllOQX, 
UK. 

A. 2 Underground facilities in USA 

• Soudan 

Latitude : 47.82° N & Longitude : 92.24° W 

Soudan Underground Lab, 30 1st Avenue, Soudan, MN 55782. 

• WIPP 

Latitude : 32.37° N & Longitude : 104.23° W 

The WIPP Experience Exhibit, U.S. Department of Energy, 4021 National Parks 
Highway, Carlsbad, New Mexico. 

• Homestake 

Latitude : 44.35° N & Longitude : 103.77° W 

Homestake Visitor Center, 160 West Main Street, Lead, SD 57754-1362. 

• SNOLAB 

Latitude : 46.47° N & Longitude : 81.19° W 
SNOLAB, Greater Sudbury, Ontario, Canada. 

• Henderson 

Latitude : 39.77° N & Longitude : 105.86° W 

Henderson Mine and Mill, 1746 County Road 202 Empire, CO 80438. 
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• Icicle Creek 

Latitude : 47.56° N & Longitude : 120.78° W 

Bridge Creek Campground, Leavenworth, Washington 98826. 

• San Jacinto 

Latitude : 33.86° N & Longitude : 116.56° W 

Mt San Jacinto State Park, 1 Tramway Road, Palm Springs, CA 92262-1827. 

• Kimballton 

Latitude : 37.37° N & Longitude : 80.67° W 
Kimballton, VA 24136. 

A. 3 Underground facilities in Europe 

• LNGS 

Latitude : 42.37° N & Longitude : 13.44° E 

Istituto Nazionale Di Fisica Nucleare - Laboratori Nazionali Del Gran Sasso-Ufficio 
Amministrativo, Strada Statale 17 Bis, L'Aquila, AQ 67100, Italy. 

• Pyhasalmi 

Latitude : 63.68° N & Longitude : 25.98° E 

86800 Pyhajarvi municipality in the south of Oulu province, Finland. 

• Slanic 

Latitude : 45.27° N & Longitude : 25.95° E 

Largest Salt mine in Europe, Prahova, Slanic, Romania. 

• Boulby 

Latitude : 54.56° N & Longitude : 0.81° W 

Boulby Potash Mine located just southeast of the village of Boulby, on the northeast 
coast of the North Yorkshire Moors in Redcar and Cleveland, England. 

• Canfranc 

Latitude : 42.76° N & Longitude : 0.51° W 

Laboratorio Subterraneo de Canfranc lies physically between New Road tunnel and 
Old Railway tunnel of Canfranc, Spain. 

• Frejus 

Latitude : 45.20° N & Longitude : 6.67° E 

Laboratoire souterrain de Modane, Carre Sciences, 73500 Modane, France. 
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• SUNLAB 

Latitude : 51.22° N & Longitude : 16.16° E 
Polkowice-Sieroszowice mine near Wroclaw in Poland. 

• Umbria 

Latitude : 42.98° N & Longitude : 12.64° E 
Umbria, IT in Italy. 

• Gran Canaria 

Latitude : 28.39° N & Longitude : 16.59° W 
Gran Canaria in Spain. 



A. 4 Underground facilities in Asia 

• CJPL 

Latitude : 28.15° N & Longitude : 101.71° E 

China JinPing Deep Underground Laboratory at Sichuan Province in China close to 
Jinping mountain. 

• Kamioka 

Latitude : 36.14° N & Longitude : 137.24° E 

Kamioka is located underground in the Mozumi Mine of the Kamioka Mining and 
Smelting Co. near the Kamioka section of the city of Hida in Gifu Prefecture, Japan. 

• Yang Yang 

Latitude : 37.77° N & Longitude : 128.89° E 

Yang Yang underground laboratory (Y2L) is located at a depth of 700 m under an 
earth overburd en in South Korea. 

• INO 

Latitude : 9.92° N & Longitude : 78.12° E 

Bodi West Hills Reserved Forest in Theni district of Tamil Nadu, India. 



B Details for all considered two-baseline combinations 



Here we first of all show the figures similar to Fig. 11 for the different accelerator laboratories 



for the sake of completeness: Fig. [T3| Fig. [T4J and Fig. [T5j Although there are quantitative 

Here the FNAL- 



differences, there are no qualitatively new insights, apart from Fig. 13 



Homestake option is shown separately, which indeed peaks at even lower 
large sin^ 26*13. 



5 GeV for 
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CPV@FNAL 
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[GeV] [GeV] 

Figure 13: CPV discovery reach at FNAL as a function of for all FNAL-options in Table [s] at the 
3ct CL for different values of (true) sin^ 2^13 (as given in the panels). Here we assume SF=1.0 and two 
50 kt detectors. We group the different baseline combinations according to the shorter baseline Li, as 
shown in the legend. The bands basically reflect the variation of the second baseline. The IDS-NF baseline 
combination is shown by the black curves with Li = 4000 km and L2 — 7500 km. 
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CPV@CERN 
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Fi glire 14: CPV discovery reach at CERN as a function of -E^i for all CERN-options in Table [s] at the 
3a CL for different values of (true) sin^ 2^13 (as given in the panels). Here we assume SF=1.0 and two 
50 kt detectors. We group the different baseline combinations according to the shorter baseline Li, as 
shown in the legend. The bands basically reflect the variation of the second baseline. The IDS-NF baseline 
combination is shown by the black curves with Li = 4000 km and L2 — 7500 km. 
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CPV@J-PARC 
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[GeV] [GeV] 

Fi glire 15: CPV discovery reach at J- PARC as a function of i?^ for all J-PARC-options in Table [s] at 
the 3cr CL for different values of (true) sin^ 26113 (as given in the panels). Here we assume SF=1.0 and two 
50 kt detectors. The IDS-NF baseline combination is shown by the black curves with Li = 4000 km and 
L2 = 7500 km. 
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Figure 16: Discovery reach for CPV, MH, and 0i3 for the combinations listed in Table [s] Dotted (red) 
curves: SF=1 (racetrack-shaped storage rings), dashed (blue) curves: SF from Table |3] (t riangular ring), 
black curves: IDS-NF baseline combination. Two 50 kt detectors used, 3ct CL. 
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Figure 17: Discovery reacli in sin^ 26'i3 as a function of baseline for CPV, MH, and ^13 



discovery and specific fractions of 6cp (contours) for a TASD, see Ref. 48 for details. The 
upper row is for one polarity stored), the lower for two polarities. The results are shown 
at 3(7 CL with true normal hierarchy. 

In Fig. [Tg] (see also following pages), we show the performances of the discussed combina- 
tions (7^1 to #50) from Table [s] In each case, we compare the performances with racetracks 
(SF=1, dotted curves), triangular-shaped geometry (SF from Table |3| dashed curves), and 
IDS combination (SF=1, solid curves). The main results of these figures are already dis- 



cussed in Sec. 5 The data points for the individual curves can be obtained at Ref. 40 



C Note on the "bi-magic" baseline 

In Ref. |48], the effect that the dependence on 5cp at a particular basehne and energy 
disappears ( "bi-magic baseline" , L ~ 2540 km) for a chosen mass hierarchy, has been studied 
for a low energy neutrino factory together with a magnetized totally active scintillator 
detector (TASD), which has somewhat different characteristics than the detector in this 
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study (somewhat better efficiencies at low energies and a better energy resolution). See 
also Ref. |51j for the bi-magic baseline in the context of a superbeam, and Ref. |52j for its 



baseline optimization. We have reproduced the simulation in Ref. 48 using the same beam 
and detector parameters, and we have studied the baseline dependence, see. Fig. 17, for the 
CPV, MH, and ^13 discovery reaches. In this figure, specific values for the fraction of 6cp 
(contours) have been chosen to show the effect of the bi-magic baseline. 



Our main observations from Fig. 17 with respect to Ref. 48 can be summarized as follows: 



1. The optimal baseline choice depends on the performance indicator, the bi-magic base- 
line is mostly preferred for the MH. Measuring CPV would clearly lead to a different 
baseline optimization. 

2. There is clear preference for using both muon polarities, which will lead to a much 
better absolute sensitivity. 

3. There is no particular preference for exactly this baseline value, for none of the per- 



formance indicators, in the sense of the "magic baseline" 11 where a clear, narrow 
dip can be seen in the numerical study (see Figs. 3 and 6 in Ref. [iS]). In addition, 
the "bi-magic effect" cannot be clearly attributed to the particular suggested energy 
windows. For example, if one masks the bins around the bi-magic energies (0.3 GeV 
windows around 1.9 and 3.3 GeV), the sensitivity is hardly affected. 

We have tested that our observations do not rely on the true hierarchy, or the energy 
resolution of the detector. 

In conclusion, the baseline of 2540 km is a good choice for the TASD if one wants to measure 
the MH as primary performance indicator for as small as possible values of sin^ 2^13. There 
is, however, no preference of this exact baseline value, instead a relatively wide baseline 
window is allowed. In addition, if CPV is considered as most important, a different baseline 
optimization will be clearly preferred. 

Finally, note that the logic of our paper is different: we have re-established the one baseline 
option of the low energy Neutrino Factory for sin^ 2^13 > 0.01. In this case, the mass 
hierarchy can be measured for all values of 6cp for L > 1 500 km for the considered TASD 
(in the most pessimistic case sin^26'i3 = 0.01). This means that the baseline optimization 
will be driven by the CPV measurement under the boundary condition of a minimal baseline 



for the mass hierarchy, which is illustrated in Figs. 1 and 3 of Ref. 37 . From Fig. 17 (lower 
left panel), we can read off that for sin^2^i3 = 0.01, the smallest considered value for 
the LENF, the optimal baseline for both MH and CPV is then 1 500 km < L < 2 100 km 
(requiring a fraction of Sqp of 80% for CPV and 100% for MH), which is roughly consistent 
to the MIND detector, see Fig. [s] upper right panel. On the other hand, for sin^ 26*13 ^ 0.01 
(^13 not discovered by the next generation of experiments), a two baseline HENF will be 
the optimal choice, and the MH discovery will be driven by the very long baseline. 
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